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Abstract 


In  this  report,  we  undertake  the  task  of  verifying  the  VLSTRACK  model  on  the  specific 
example  of  raindrop  terminal  velocities,  which  were  calculated  directly  from  Fluid  Dynamics. 
In  this  example,  the  verification  is  in  good  agreement  with  experiments. 
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Executive  Summary 


In  this  research,  we  undertook  a  task  to  verify  the  validity  of  the  Vapor,  Liquid,  and  Solid 
Tracking  (VLSTRACK)  Computer  Model  in  calculating  terminal  velocity  of  raindrops.  The  findings 
compare  very  well  with  Fluid  Dynamics  calculations  and  experiments. 

The  version  VLSTRACK  1.6.2,  used  here,  is  developed  and  owned  by  the  U.S.  Naval  Surface 
Warfare  Center  (NSWC),  Dahlgren  Division,  Dahlgren,  VA  22448-5100. 

If  the  revision  of  VLSTRACK  is  contemplated,  we  recommend  that  more  detailed  documentation 
of  input  parameters  for  explaining  specific  tasks  be  included. 


Vll 


Intentionally  left  blank. 


1.  Introduction 


In  this  report,  we  are  verifying  whether  the  VLSTRACK  [1]  computer  model  is  applicable  for 
a  specific  example — comparison  of  terminal  velocities  of  falling  water  droplets  under  the  influence 
of  gravity  and  drag  forces  (the  drag  forces  contain  ordinary  drag  forces  and  buoyancy  forces). 
Specifically,  we  are  comparing  the  VLSTRACK  calculated  terminal  velocities  in  the  atmosphere 
with  those  calculated  interactively  through  fluid  dynamics  formalism  [2]. 

Although  one  might  think  that  this  is  a  simple  task,  it  turned  out  that  before  one  could  set 
VLSTRACK  to  calculate  the  terminal  velocities,  a  lot  of  preliminary  manipulations  were  necessary 
in  the  data  input  files  in  order  to  successfully  accomplish  this  comparison. 

Section  2  outlines  the  purpose  of  this  analysis.  Section  3  is  devoted  to  recapitulation  of  fluid 
dynamics  calculations  of  raindrop  terminal  velocities  and  the  experimentally  determined  values.  The 
analysis  using  VLSTRACK  is  discussed  in  section  4.  Section  5  provides  the  velocities  and 
summarizes  the  results.  The  significant  results  are  discussed  in  section  6. 

The  purpose  of  this  report  is  to  show  that  an  accurate  description  of  a  particular  situation  by 
VLSTRACK  requires  significant  manipulation  of  input  VLSTRACK  parameters.  Our  specific 
example  of  VLSTRACK  calculated  raindrop  terminal  velocities  and  their  subsequent  comparison 
with  fluid  dynamics  calculations  [2]  and  experimental  data  [3]  illustrates  how  comparisons  for  other 
agents  may  be  made. 

2.  Theory:  Overview  of  Fluid  Dynamics  Calculations  of 
Raindrop  Terminal  Velocities 

Here  we  give  a  short  review  of  fluid  dynamics  calculations  of  raindrop  terminal  velocities.  A 
full  account  may  be  found  in  Soln  [2]. 
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When  a  raindrop  is  of  a  diameter  that  is  smaller  than  5  mm,  one  can  treat  a  freely  falling  droplet 
as  a  sphere;  a  fact  that  can  be  verified  a  posteriori  by  comparison  with  experiments. 

Treating  a  droplet  as  a  sphere  with  water  density  simplifies  the  calculation  of  the  terminal 
velocities. 

For  the  sake  of  completeness,  we  write  down  the  differential  equation  for  the  water  sphere  when 
forces  acting  on  it  are  all  in  the  vertical  direction.  The  velocity  of  the  sphere  in  the  horizontal 
direction  is  assumed  to  be  zero.  Hence  the  differential  equation  for  the  vertical  motion  of  the  sphere 
reads  as  follows: 


d2y 

dt2 


fl-f' 

-  (1  -  P)g-  l£cd  u |u | 

\  2 ) 
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where  u  is  the  velocity  of  the  sphere  in  the  vertical  direction,  whose  coordinate  here  is  denoted  as 
y,  looking  upward. 

The  relevant  parameters  appearing  in  equation  (1)  are: 

p  =  1000.0  kg/m3,  p,  =  1.22  kg/m3,  p  =  -^  =  0.0012,  g  =  9.8  m/s2,  (2) 

P 

where  p  is  the  density  of  water,  pf  is  the  density  of  air,  and  cd  is  a  drag  coefficient  that  depends  on 
the  dimensionless  Reynolds  number 

c<j  =  cd  (Re) .  (3) 


The  Reynolds  number  has  to  be  computed  separately  for  each  velocity,  u,  according  to 
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where  d  is  the  diameter  of  the  sphere  and  v  is  the  ordinary  viscosity  of  air: 


v  =  0.0000149  m2/s . 

For  air  as  a  medium,  empirical  relation  (3)  can  be  specified  as  [2, 3] 

cH  =  1  <  R.  £  400. 

d  r,  0.646  e 


(5) 


(6) 


The  terminal  velocity  ut  follows  when  equation  (1),  acceleration,  is  equivalent  to  zero,  i.e.,  when 
the  velocity  in  the  downward  direction  becomes  constant.  Thus,  equating  (1)  to  zero  and  setting 
u  =  -u,,  where  ut  looks  in  the  downward  direction,  we  obtain: 


ut  = 


N 


4g(l-p)d 


Cd  3p 


(7) 


The  next  thing  one  does  is  to  “guess”  on  the  initial  value  of  terminal  velocity  utfor  a  droplet  with 
a  given  diameter  d.  For  example,  for  d  =  0.001  m,  we  guess  the  initial  value  for  ut  to  be  3  m/s. 
Putting  this  into  relation  (4)  for  the  velocity  u,  we  obtain  Re  =  201 .34.  This  now  is  used  to  calculate 
cd  according  to  relation  (6),  resulting  in  cd  =  0.78.  Now  we  have  everything  for  calculating  the 
“improved”  value  for  ut  according  to  relation  (7).  The  improved  value  is  ut  =  3.73  m/s.  Now  one 
repeats  the  cycle  with  this  new  value  for  ut  yielding  even  better  approximation  for  ut.  Eventually, 
this  repetition  will  yield  ut’s  that  do  not  differ  from  each  other,  which  means  that  we  found  the  most 
accurate  value  for  ut.  In  our  case,  ut  =  4.1  m/s  as  a  final  value. 

One  can  do  this  type  of  calculation  of  ut  for  every  diameter  d.  The  results  are  exhibited  and 
compared  to  equivalent  values  and  VLSTRACK  calculations  in  Table  1. 
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Table  1.  Comparison  of  VLSTRACK  Calculated,  Fluid  Dynamics  Calculated, 
and  Experimental  Terminal  Velocities  of  Water  Droplets 


Mass  Median 
Diameter  (MMD) 

(P) 

Fluid  Dynamics 
Calculations  of 
Terminal  Velocity  of 
Water  (Sphere) 
(m/min) 

VLSTRACK 
Calculations  of 
Terminal  Velocity  of 
Water  (Sphere) 
(m/min) 

Measured  Terminal 
Velocity  of  Water 
(Droplet) 
(m/min) 

500 

108.0 

122.8 

120.0 

1,000 

246.0 

232.9 

240.0 

1,500 

336.0 

305.6 

330.0 

2,000 

390.0 

379.8 

396.0 

2,500 

444.0 

407.5 

450.0 

3,000 

480.0 

431.6 

492.0 

3,500 

522.0 

448.2 

522.0 

4,000 

558.0 

453.8 

540.0 

4,500 

588.0 

458.8 

546.0 

5,000 

624.0 

463.3 

546.0 

5,500  | 

648.0 

467.4 

552.0 

6,000 

684.0 

471.2 

552.0 

The  measured  terminal  velocities  of  water  droplets  come  from  Blanchard  [3] .  The  VLSTRACK 
data  will  be  discussed  in  the  next  section. 

3.  Implementation:  Analysis  Using  VLSTRACK 

Now  that  we  have  provided  the  rationale  for  calculated  terminal  velocity,  we  want  to  look  at 
VLSTRACK  computer  model  [1]  calculated  terminal  velocities  compared  to  those  we  have 
calculated  interactively  through  the  fluid  dynamics  equation  sequence  discussed  in  Section  2. 

The  sequence  we  used  to  set  up  the  VLSTRACK  files  is  discussed  in  this  section.  We  are  going 
to  specify  the  process  used  to  compute  the  VLSTRACK  results  in  Table  1.  Before  starting  to 


interrogate  VLS TRACK,  we  must  determine  the  boundaries  of  what  comparisons  we  want  to  make. 
We  want  to  compare  VLSTRACK  calculated  values  to  the  fluid  dynamics  calculated  values,  so  we 
need  to  duplicate  the  parameters  used  in  the  fluid  dynamics  terminal  velocity  calculations  as  closely 
as  possible.  This  means  that  we  want  to  look  at  the  water  droplet  as  a  rigid  sphere.  We  therefore 
set  forth  to  define  a  VLSTRACK  “threat”  scenario  that  will  replicate  the  physical  conditions  and 
properties  of  a  rigid  sphere  of  water  and  will  mimic  stabile  atmospheric  conditions. 

The  first  step  in  this  definition  sequence  was  to  establish  an  agent  in  the  agent  parameter  file 
named  VLSAGN.PAR  [1]  to  define  the  physical  properties  of  the  rigid  sphere  of  water.  The  agent 
file  was  modified  to  identify  water  as  a  biological  type  agent  (type  5),  meaning  that  the  water  would 
be  treated  as  particle  beads  rather  than  as  a  liquid. 

Agent  parameter  file  specifications  also  include  the  bulk  density,  dissemination  efficiency, 
median  lethal  or  effective  dosage,  probit  slope,  and  freezing  temperature.  Dissemination  efficiency 
was  set  at  100%,  and  median  lethal  or  effective  dosage  was  set  at  only  0.001  mg-mi/m3  to  assure  the 
water  would  exhibit  deposition  (in  the  case  of  deposition  runs)  and  that  effects  would  be  detected 
(in  the  case  of  dosage  effects  runs).  The  bulk  density,  probit  slope,  and  freezing  temperature  were 
set  equivalent  to  that  of  water.  In  addition  to  the  overall  agent  identification,  each  type  5  agent 
includes  particle  MMD  [1],  geometric  droplet/particle  distribution  sigma,  particle  density,  daytime 
biological  decay  rate,  nighttime  biological  decay  rate,  biological  agent  purity,  and  biological  agent 
dissemination  efficiency  for  exploding  munitions.  The  particle  MMD  was  varied  from  500  p  to 
6,000  p  in  increments  of  500  p .  The  fluid  dynamics  calculations  are  valid  up  to  a  water  droplet 
diameter  of 5,000  p .  We  included  up  to  6,000  p  diameters  in  order  to  exhibit  explicitly  the  departure 
from  experimental  values.  The  geometric  droplet/particle  distribution  sigma  was  set  at  1,  signifying 
that  all  droplets  are  the  same  size;  the  particle  density  was  set  at  that  of  water  (at  21  °  C)  equal  to  1 
g/cm3;  and  the  biological  agent  dissemination  efficiency  is  not  pertinent  since  it  is  not  used  when 
user-defined  (i.e.,  puff-type,  localized  concentration  of  a  gaseous  or  liquid  substance  [  1  ])  munitions 
are  utilized.  Daytime  and  nighttime  biological  decay  rates  were  set  at  zero,  and  the  biological  agent 
purity  was  set  at  100%  to  maximize  capacity  of  the  water  “hazard.” 
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The  second  sequential  step  was  to  define  a  munition  in  the  parameter  file  named  VLSMUN.PAR 
[1],  which  is  the  location  where  munition  descriptions  are  built  or  modified.  It  was  decided  to 
disseminate  water  in  puffs  as  a  line  source  to  view  the  water  from  each  puff  as  it  was  deposited  via 
the  graphics  output  file,  to  track  and  compare  terminal  velocities  of  each  puff  by  reviewing  the 
output  files  named  VLSTRACK.TRK  and  VLSTRACK.REC,  respectively  [1].  The  munition  was 
categorized  as  a  user  defined  munition  and  then  further  defined  in  a  separate  puff  property  file  named 
VLSTRACK.POS  [1].  This  file  was  configured  to  contain  five  records  for  each  munition  puff, 
including  the  mass,  downrange  and  crossrange  distance,  height,  and  detonation  time.  After  making 
several  successful  VLSTRACK  runs,  in  terms  of  measuring  terminal  velocity,  a  standard  scenario 
was  established  and  then  the  puff  property  file  parameters  were  varied  many  times.  Keeping  in  mind 
that  it  was  desired  to  make  the  puff  release  height  higher  than  the  distance  the  droplet  would  travel 
in  1  min,  the  puff  parameter  values  were  varied  to  substantiate  uniformity  of  terminal  velocity  from 
run  to  run  and  to  determine  terminal  velocity  for  each  different  MMD. 

The  third  step  was  to  develop  a  meteorology  file  in  a  file  named  VLSTRACK.MET.  We  know 
that  meteorological  conditions  have  a  significant  impact  on  how  agents  transport  and  disperse  on  the 
battlefield.  The  object  of  our  VLSTRACK.MET  file  is  to  define  atmospheric  conditions  that  parallel 
the  constant  conditions  used  for  the  fluid  dynamics  calculations  [2]  for  the  droplets  under  the 
influence  of  only  gravity  and  drag.  Preliminary  VLSTRACK  runs  were  performed  using  a  time 
series  meteorology  file  in  which  meteorology  data  was  input  for  specific  time  periods  at  one 
measured  height.  The  model  computes  meteorological  data  for  all  other  pertinent  heights  based  on 
the  meteorological  input  data  and  on  the  atmospheric  stability  category  known  as  the  Pasquill 
category  [1].  VLSTRACK  is  programmed  to  either  allow  the  program  to  specify  the  Pasquill 
stability  category  based  on  the  meteorological  conditions,  time  of  day,  and  terrain  information  or  to 
allow  the  operator  to  select  a  particular  stability  category.  There  are  seven  categories  available, 
ranging  from  very  unstable  to  very  stable.  In  this  case,  the  option  was  selected  for  the  program  to 
determine  the  Pasquill  category. 

A  height-time  variable  meteorology  file  was  used  for  all  MMD  comparative  analysis  runs  in 
order  to  fix  the  wind  at  a  constant  0  km/hr  and  the  temperature  at  a  constant  21  °  C  at  all  elevations 
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to  parallel  the  conditions  of  the  fluid  dynamics  calculations  [2].  It  was  necessary  to  use  the  0.001- 
km/hr  wind  speed  because  an  anomaly  of  the  model  did  not  allow  a  constant  0  km/hr  wind  speed. 
The  height-time  variable  meteorology  file  includes  ground  surface  type,  Pasquill  stability  category, 
cloud  cover,  terrain  elevation  above  sea  level,  threat  measurement  height  above  sea  level,  wind 
bearing,  wind  speed,  air  temperature,  and  atmospheric  pressure.  The  height-time  variable 
meteorology  file  that  was  used  consisted  of  two  identical  time  steps.  Including  the  constant 
temperature  and  wind,  the  additional  parameters  were  brush-type  ground  surface,  neutral  Pasquill 
stability  category,  clear — no  cloud  cover,  and  0-m  terrain  elevation  above  sea  level  for  location  to 
look  at  the  deposition.  Four  records  were  included  for  each  time  period.  They  were  as  follows: 

•  0.001 -km  measurement  height  above  sea  level,  90°  true  north  (DTN)  [1]  wind  bearing, 
0.001-km/hr  wind  speed,  21 0  C  air  temperature  at  1,000-mbar  atmospheric  pressure. 

•  0.5-km  measurement  height  above  sea  level,  90  DTN  wind  bearing,  0.001-km/hr  wind  speed, 
21  °  C  air  temperature  at  950-mbar  atmospheric  pressure. 

•  1 .0-km  measurement  height  above  sea  level,  90  DTN  wind  bearing,  0.001-km/hr  wind  speed, 
21  °  C  air  temperature  at  900-mbar  atmospheric  pressure. 

•  2.0-km  measurement  height  above  sea  level,  90  DTN  wind  bearing,  0.001-km/hr  wind  speed, 
21 0  C  air  temperature  at  800-mbar  atmospheric  pressure. 

After  the  water  agent,  puff  munitions,  and  meteorological  data  were  detailed  in  the  appropriate 
parameter  files,  the  next  step  was  to  complete  an  attack  scenario  using  the  six  available  VLSTRACK 
editor  windows.  Included  is  the  time  of  day  for  onset  of  the  attack — in  this  case,  release  of  the  first 
puff.  All  of  the  modeling  scenarios  were  performed  using  0400  (military  denoted)  local  onset  time 
since  night  time  and  early  morning  conditions  provide  the  calmest  atmospheric  conditions  [4]. 
Cumulative  deposition  intervals  of  1  min  were  used  for  run  output.  Level  of  detail  of  deposition 
contour  levels  was  set  as  low  as  possible — 0.01  mg/m2  in  order  to  detect  any  deposition. 
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4.  Example  Application:  Number  Crunching  of  Terminal 
Velocity 

Model  scenarios  were  built  and  run  one  at  a  time,  with  each  subsequent  run  description 
contingent  on  the  preceding  run  results.  Output  graphics  of  hazard  footprints  were  viewed  at  each 
interval  during  each  model  run  to  monitor  the  time  at  which  water  deposited  on  the  ground.  Model 
scenario  inputs  were  modified  until  run  output  graphics  exhibited  deposition.  A  uniform  base 
scenario  was  then  defined  and  mn  sequentially  varying  only  sphere  MMDs  from  500  p  to  6,000  p. 
Terminal  velocity  numbers  were  calculated  by  looking  at  each  individual  munition  puff  and  by 
analyzing  all  resulting  output  files.  In  particular,  after  viewing  deposition  contours  in  the  graphics 
output  display  file,  the  corresponding  VLSTRACK.TRK  output  file  was  saved  for  each  file.  This 
VLSTRACK.TRK  file  shows  tracking  for  each  droplet  (discharged  in  the  munition  puff),  including 
such  data  as  the  x  and  y  distance  traveled  and  the  height  of  the  droplet  at  that  record.  The  terminal 
velocity  was  approximated  by  looking  at  the  time  it  took  for  deposition  to  start  and  the  height  the 
droplets  were  released  as  recorded  in  the  output  file  VLSTRACK.OUT  file  [1].  The 
VLSTRACK.TRK  file  was  also  used  to  calculate  an  exact  terminal  velocity  by  calculating  the 
difference  in  the  height  of  the  droplet  from  one  time  snapshot  to  the  next. 

5.  Conclusions  and  Results 

The  VLSTRACK.TRK,  Table  2,  and  the  VLSTRACK.REC,  Table  3,  output  files  track  and 
record  the  cloud  properties  of  all  airborne  vapor,  droplet,  or  particle  clouds  for  each  individual  puff 
released.  These  files  were  used  to  calculate  and  compare  the  terminal  velocities  of  water  beads  from 
each  puff.  The  resulting  terminal  velocity  calculations  from  tracked  puffs  were  all  consistent  for 
each  individual  MMD. 

The  results  of  the  VLSTRACK  calculated  terminal  velocities  were  within  expected  tolerance  of 
the  fluid  dynamics  calculations  up  to  a  MMD  of 2,500  p .  Table  1  shows  the  VLSTRACK  and  fluid 
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dynamics  [2]  calculated  terminal  velocity  of  a  rigid  sphere  of  water  and  the  measured  terminal 
velocity  of  rain  droplets. 


The  fluid  dynamics  calculations  are  consistent  with  the  measured  terminal  velocities  of  up  to 
5,000  p  MMD,  while  VLSTRACK  calculations  are  consistent  to  water  droplets  of  up  to  2,500  p 
MMD.  In  fact  this  conclusion  can  also  be  obtained  through  statistical  analysis  which  compares 
calculated  values  with  each  other  and  with  experiments  [5]. 


6.  Discussion  of  the  Significant  Results 


The  final  results  of  this  study  showed  that  VLSTRACK  models  terminal  velocity  accurately  for 
rigid  spheres  of  water  with  MMDs  from  500  to  2,500  p  as  indicated  by  the  heavy  line  in  Table  1. 
The  larger  the  MMD,  >2,500  p,  the  more  discrepancy  in  the  results.  These  results  are  reassuring 
and  could  be  expected  since  chemical  agents  that  VLSTRACK  is  intended  to  model  would  not  have 
such  large  MMDs  anyway.  For  example,  typical  MMDs  for  most  dusty  and  bio  agents  are  only  5  p , 
for  glass  beads  are  250  p,  for  neat  agent  droplets  are  500  p,  and  for  thickened  agent  droplets  are 
2,500  p  maximum.  VLSTRACK  terminal  velocity  calculations  also  remained  uniform  for  each 
different  MMD  even  when  altering  mass  or  release  position.  The  discrepancy  of  the  fluid  dynamics 
calculations  to  the  measured  terminal  velocities  beyond  5,000  p  MMD  would  indicate  a  breakup  of 
droplets. 

This  comparison  study  clearly  indicates  the  validity  of  VLSTRACK  for  evaluating  terminal 
velocities  in  the  range  for  which  it  is  intended. 
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Table  2.  Example  of  the  VLSTRACK.TRK  Output  File  Containing  the  Properties  of  the  Droplets  as  They  Propagate  Through 
the  Atmosphere 
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Note:  The  columns  are  for  time  in  seconds,  droplet  size  group  number,  puff  number,  downwind  distance  in  meters,  crosswind  distance  in  meters,  height  above  sea  level  in  meters, 
downwind  cloud  sigma  in  meters,  crosswind  cloud  sigma  in  meters,  vertical  cloud  sigma  in  meters,  average  terrain  elevation  in  meters,  local  terrain  elevation  in  meters,  mixing  layer 
height  in  meters,  downwind  cloud  rotation  angle  in  degrees,  droplet  radius  in  microns,  cloud  concentration  in  milligrams  per  cubic  meter,  and  cloud  mass  in  milligrams,  respectively. 


Table  2.  Example  of  the  VLSTRACK.TRK  Output  File  Containing  the  Properties  of  the  Droplets  as  They  Propagate  Through 
the  Atmosphere  (continued) 


Note:  The  columns  are  for  time  in  seconds,  droplet  size  group  number,  puff  number,  downwind  distance  in  meters,  crosswind  distance  in  meters,  height  above  sea  level  in  meters, 
downwind  cloud  sigma  in  meters,  crosswind  cloud  sigma  in  meters,  vertical  cloud  sigma  in  meters,  average  terrain  elevation  in  meters,  local  terrain  elevation  in  meters,  mixing  layer 
height  in  meters,  downwind  cloud  rotation  angle  in  degrees,  droplet  radius  in  microns,  cloud  concentration  in  milligrams  per  cubic  meter,  and  cloud  mass  in  milligrams,  respectively. 


Table  2.  Example  of  the  VLSTRACK.TRK  Output  File  Containing  the  Properties  of  the  Droplets  as  They  Propagate  Through 
the  Atmosphere  (continued) 


Note:  The  columns  are  for  time  in  seconds,  droplet  size  group  number,  puff  number,  downwind  distance  in  meters,  crosswind  distance  in  meters,  height  above  sea  level  in  meters, 
downwind  cloud  sigma  in  meters,  crosswind  cloud  sigma  in  meters,  vertical  cloud  sigma  in  meters,  average  terrain  elevation  in  meters,  local  terrain  elevation  in  meters,  mixing  layer 
height  in  meters,  downwind  cloud  rotation  angle  in  degrees,  droplet  radius  in  microns,  cloud  concentration  in  milligrams  per  cubic  meter,  and  cloud  mass  in  milligrams,  respectively. 


Table  2.  Example  of  the  VLSTRACK.TRK  Output  File  Containing  the  Properties  of  the  Droplets  as  They  Propagate  Through 
the  Atmosphere  (continued) 
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Note:  The  columns  are  for  time  in  seconds,  droplet  size  group  number,  puff  number,  downwind  distance  in  meters,  crosswind  distance  in  meters,  height  above  sea  level  in  meters, 
downwind  cloud  sigma  in  meters,  crosswind  cloud  sigma  in  meters,  vertical  cloud  sigma  in  meters,  average  terrain  elevation  in  meters,  local  terrain  elevation  in  meters,  mixing  layer 
height  in  meters,  downwind  cloud  rotation  angle  in  degrees,  droplet  radius  in  microns,  cloud  concentration  in  milligrams  per  cubic  meter,  and  cloud  mass  in  milligrams,  respectively. 


Table  2.  Example  of  the  VLSTRACK.TRK  Output  File  Containing  the  Properties  of  the  Droplets  as  They  Propagate  Through 
the  Atmosphere  (continued) 


14 


Note:  'Hie  columns  are  for  time  in  seconds,  droplet  size  group  number,  puff  number,  downwind  distance  in  meters,  crosswind  distance  in  meters,  height  above  sea  level  in  meters, 
downwind  cloud  sigma  in  meters,  crosswind  cloud  sigma  in  meters,  vertical  cloud  sigma  in  meters,  average  terrain  elevation  in  meters,  local  terrain  elevation  in  meters,  mixing  layer 
height  in  meters,  downwind  cloud  rotation  angle  in  degrees,  droplet  radius  in  microns,  cloud  concentration  in  milligrams  per  cubic  meter,  and  cloud  mass  in  milligrams,  respectively. 


Table  2.  Example  of  the  VLSTRACK.TRK  Output  File  Containing  the  Properties  of  the  Droplets  as  They  Propagate  Through 
the  Atmosphere  (continued) 


15 


Table  2.  Example  of  the  VLSTRACK.TRK  Output  File  Containing  the  Properties  of  the  Droplets  as  They  Propagate  Through 
the  Atmosphere  (continued) 
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Note:  The  columns  are  for  time  in  seconds,  droplet  size  group  number,  puff  number,  downwind  distance  in  meters,  crosswind  distance  in  meters,  height  above  sea  level  in  meters, 
downwind  cloud  sigma  in  meters,  crosswind  cloud  sigma  in  meters,  vertical  cloud  sigma  in  meters,  average  terrain  elevation  in  meters,  local  terrain  elevation  in  meters,  mixing  layer 
height  in  meters,  downwind  cloud  rotation  angle  in  degrees,  droplet  radius  in  microns,  cloud  concentration  in  milligrams  per  cubic  meter,  and  cloud  mass  in  milligrams,  respectively. 


Table  2.  Example  of  the  VLSTRACK.TRK  Output  File  Containing  the  Properties  of  the  Droplets  as  They  Propagate  Through 
the  Atmosphere  (continued) 
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Note:  The  columns  are  for  time  in  seconds,  droplet  size  group  number,  puff  number,  downwind  distance  in  meters,  crosswind  distance  in  meters,  height  above  sea  level  in  metere, 
downwind  cloud  sigma  in  meters,  crosswind  cloud  sigma  in  meters,  vertical  cloud  sigma  in  meters,  average  terrain  elevation  in  meters,  local  terrain  elevation  in  meters,  mixing  layer 
height  in  meters,  downwind  cloud  rotation  angle  in  degrees,  droplet  radius  in  microns,  cloud  concentration  in  milligrams  per  cubic  meter,  and  cloud  mass  in  milligrams,  respectively. 


Table  2.  Example  of  the  VLSTRACK.TRK  Output  File  Containing  the  Properties  of  the  Droplets  as  They  Propagate  Through 
the  Atmosphere  (continued) 
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Note:  The  columns  are  for  time  in  seconds,  droplet  size  group  number,  puff  number,  downwind  distance  in  meters,  crosswind  distance  in  meters,  height  above  sea  level  in  meters, 
downwind  cloud  sigma  in  meters,  crosswind  cloud  sigma  in  meters,  vertical  cloud  sigma  in  meters,  average  terrain  elevation  in  meters,  local  terrain  elevation  in  meters,  mixing  layer 
height  in  meters,  downwind  cloud  rotation  angle  in  degrees,  droplet  radius  in  microns,  cloud  concentration  in  milligrams  per  cubic  meter,  and  cloud  mass  in  milligrams,  respectively. 


Table  3.  Example  of  the  VLSTRACK.REC  Output  File  Containing  the  Properties  of  the  Droplets  Impacting  the  Ground  Surface 
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Note:  The  columns  are  for  cloud  smoothing  index,  droplet  size  group  number,  puff  number,  time  in  seconds,  droplet  diameter  in  microns,  downwind  distance  in 
meters,  cross  wind  distance  in  meters,  downwind  cloud  sigma  in  meters,  cross  wind  cloud  sigma  in  meters,  downwind  cloud  rotation  angle  in  degrees,  and  concentration 
in  milligrams  per  cubic  meter,  respectively. 


Table  3.  Example  of  the  VLSTRACK.REC  Output  File  Containing  the  Properties  of  the  Droplets  Impacting  the  Ground 
Surface  (continued) 
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